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Abstract 

Assume (1) asset returns follow a stochastic multi-factor process with 
time-varying conditional expectations; (2) investments are linear func- 
tions of factors. This paper calculates asymptotic joint moments of the 
logarithm of investor's wealth and the factors. These formulas enable fast 
computation of a wide range of investment criteria. The results are illus- 
trated by a numerical example that shows that the optimal portfolio rules 
are sensitive to the specification of the investment criterion. 

1 Introduction 



The asset returns are predictable: see (Cochrane 1999) and (Campbell, Lo 
& MacKinlay 1997) who overview the empirical support for asset return pre- 
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dictability^, and (Brennan, Schwartz & Lagnado 1997), (Campbell & Viceira 
1999), and (Bielecki & Pliska 1999) who study optimal portfolio allocation in 
the situation with predictable returns. Conventional dynamic stochastic pro- 
gramming needs a lot of computing power to find optimal allocations, and a 
possible approach to this problem is to choose an approximate investment cri- 
terion that would make the problem numerically tractable. This paper studies 
portfolio optimization for a class of such criteria and tests the sensitivity of 
optimal portfolio rules to the choice of the investment criterion. 

In early 1970s Samuelson and Merton derived Hamilton- Jacobi-Bellman equa- 
tion for the solution of the dynamic asset allocation problem, an equation that 
is valid for any type of return process and any type of investment criteria.^ 
This partial differential equation is, however, hard to solve. To circumvent this 
difficulty, (Bielecki & Pliska 1999)^ introduced the risk-sensitive investment cri- 
terion, which makes the optimization problem easier. The criterion assumes 
that the investor's preferences depend only on the growth rate of the expected 
logarithm of portfolio value and on the growth rate of its variance. This assump- 
tion appears to be excessively tight since it ignores other relevant moments, such 
as the covariance with the factors. An investor might be concerned with the 

performance of the portfolio in "bad" times, when his other sources of income 
^For example, see empirical studies by (Balvers, Cosimano & McDonald 1990), (Breen, 

Glosten & Jagannathan 1990), (Campbell 1987), (Cochrane 1991), (Fama & French 1989), 

(Pcsaran & Timmcrmann 1995), (Pcsaran & Timmcrmann 2000), (Cooper 1999). 
2 (Samuelson 1969), (Merton 1969), (Merton 1971), and (Merton 1973) 
^See also (Stettner 1999), (Bielecki, Pliska & Sherris 2000), (Bielecki, Hernandez-Hernandez 

& PUska 1999), (Bielecki & Pliska 2001) 
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are low. The criterion that depends only on variance and expectation of the 
portfolio value ignores this consideration. 

This paper uses a new method to compute the asymptotic joint moments of 
the logarithm of the portfolio value and factors, allowing for the fast compu- 
tation of a wide class of investment criteria. An example calibrated on the real 
data shows that the optimal portfolio rules strongly depend on the specification 
of the investment criterion. In particular, the rule is sensitive to the inclusion of 
the correlation of portfolio value with factors in the criterion. This supports the 
view that the correct choice of the approximate investment criterion is crucial 
for optimal portfolio allocation. 

The remainder is organized as follows. Section 2 explains the model. Section 
3 formulates the main result. Section 4 specializes the result to a one-factor 
example, calibrates it on the real data and presents numerical illustrations. 
And Section 5 concludes. The proof of the main result is in Appendix. 



2 Model 

Assumption 1. Securities follow a dynamic factor model. 



-^ = {a + AX{t))idt+^aikdWk{t), i=l,2,...,m, (1) 

^'''^> k=i 

dX{t) = BX{t)dt + MW{t), (2) 

where W{t) is a 3R™+" valued standard Brownian motion process, X{t) is 5R" 
valued factor process. It is assumed that nxn matrix B is stable. This coincides 



with the model of (Bielecki & PHska 1999) and (Bielecki et al. 2000) with the 
exception that the factor process X{t) is normalized to have zero mean. 
The process for investor's wealth U is 

dU = ^dS, (3) 

where / is the investment vector measured in shares of wealth. 
Assumption 2. / is a linear function of factors. 

I = h + HX, (4) 

where H is a. constant mx n matrix and /i is a constant m x 1 vector.* 
Thus, the investor's wealth U follows the process 

7 Q 

dU = U{h + HXy— (5) 

Assumption 3. The investment criterion is 

W{h, H) = liminf -{Eu{t) - ^var{u{t))} + T lim E{u{t)X{t)), (6) 

where u{t) =: ln[/(i) - ln[/(0). 

Note that this criterion includes the covariance of the logarithm of portfolio 
value with the factors. An interpretation of this criterion is that the investor 
cares not only about growth of the portfolio value and its volatility but also 
about the ability of the portfolio to generate good returns in bad times when 
his other income sources bring low returns. 



^This is not a very stringent assumption because non-linear functions of the factors can 
always be added to the set of all factors to incorporate non-linear dependence of investments 
on the factors. 



3 Result 

The following theorem shows how to compute asymptotic moments that enter 
the investment criterion. 

Theorem 1 If assumptions 1 and 2 hold then 

E{u) = t[h'a - ^h'EE'h + tr{A{H'A - ^H'EE'H))], (7) 
Var(w) ~ t[YY' + tr{2SH'A + (A - S)H'^^'H)] + const as t ^ oo, (8) 
E{uX) ~ B-^ [A(iJ'SS'/i - A'h - H'a) - AS'/i] as t ^ oo. (9) 

Here A^ =: E{xiXj), vector Y is defined as follows: 



Y =: (/I'SE'fl- - h'A - a'H)B-^K + h'T,, (10) 

and matrix S is the symmetric solution of the equation 

BS + SB' = -2AiI'^A + AH'HH'HA - 2AS'iJA. (11) 

The proof of Theorem 1 is in the Appendix. Let me explain here its main 
idea for the case with only one factor x. Suppose that we want to compute 
limt^oo £(m*x^) and we already know limt^^ E{u''x'') ior k < i and for k = 
i,l < j. Using Ito Lemma, we can write 

d{u^x^) =iu'~^x^du + jx^~^u'dx 

+ ^{*(« - l)u'-^x\duy + j{j - l)u'x^-^{dxf + iju'-^x^-^dudx). 

(12) 
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After the expressions for du and dx are substituted and the expectations are 
taken, this equation is reduced to the following differential equation on E{u^x^): 



where . . . includes only terms with already known asymptotic limits. Since B 
is stable (in this one factor example it means B <Q), the solution of this differ- 
ential equation has an asymptotic limit that can be easily calculated. Theorem 

1 extends this method to the multi-factor case. 

4 Application 

This section applies the general result to the particular example with the param- 
eters estimated from real data. First, Theorem 1 is specialized to the case with 
only one asset and one factor. Thus, a. A, B, h, and H are scalars, S = {a, ij), 
and A = (0, A). 

In this setting. Theorem 1 is reduced to the following 

Corollary 1 



—E{u'x^) = jBE{u'x^) + ... 



dt 



(13) 



E{u) = t{ha - ^HA - ^^{h' - ^H^)} 



(14) 



E{ux) ^ ^{hHiri" + a^) -Ah- Ha} - -^h 
Var{u) ~ t{YY' + 2SHA + (^ - S)H^{a^ + if)} 



(15) 



(16) 



where 



Y 



{hH{a^ + rf) -hA- aH}hQ, A) + h{cT, rj) 




S 
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I have estimated the parameters of the stochastic process using the monthly 
data from CRSP files. The effective stock return is taken as the difference of the 
rate of return on S&P composite index and the short-term interest rate. The 
factor is the short-term interest rate, which is the rate on 3-month Treasury 
Bill. The results of estimation are summarized in the next two tables. 



Table 1. Process Estimates" 





constant 


r 


dS/S 


1.993 
(3.505) 


-1.177 
(-14.220) 


dr 


0.120 
(0.911) 


0.979 
(42.885) 



''Period is from January 1970 to December 2000. Rates and returns are measured in percentage terms; t-ratios are 
in parentheses 

Table 2. Covariances of interest rate and stock return innovations" 





Slock Rcluni 


r 


Stock Return 


19.587 


0.0553 


r 


0.0553 


0.4006 



''Period is from January 1970 to December 2000 



Therefore, the parameters of the model can be calibrated as a = 0.01993, 
A = -0.01177, B = -0.021, A = 0.6329, r] = 0.000874, and a = 0.044249. 

For the first illustration, assume that /i = 1, so that on average the investor 
have all his funds in the stock. If the interest rate deviates from its average value 
then the investor can move his assets to bonds. He can also borrow additional 
funds and increase his position in the stock. The first three graphs show how 
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Eu, Var{u), and E{ux) change with respect to changes in H. 

These graphs show that the strategy of investing more when the interest rates 
are low increases the expected return of the portfoho. Expected covariance of 
portfoho value with factors increases for large H. Surprisingly, the volatility of 
the portfolio is minimized at a certain non-zero level of H. 

For the second illustration, assume that the target function (6) is specified 
with factor-sensitivity F = 0. Then the dependence of optimal [h, H) combina- 
tion on the risk-sensitivity parameter 6 is illustrated on Figure 4. 

This picture shows that the increase in the risk-sensitivity leads to smaller 
average investment in stocks and to smaller sensitivity of the investment to 
changes in the interest rates. It is interesting to note, however, that the ratio 
of H to /i, which can be interpreted as the relative sensitivity of investments to 
interest rates is almost constant for sufficiently large values of risk-sensitivity. 

Now consider factor-sensitivity parameter F changing from to 0.01. The 
graph on Figure 5 shows the optimal (/i, H) path. The greater sensitivity of 
the investment rule to the covariance of the portfolio value with the factor leads 
to smaller average investment and smaller sensitivity of the investment rule to 
the factor. It is interesting, however, that the sensitivity of investment rule to 
factor is affected smaller than the average investment amount. 
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Figure 1: Dependence of K = E{du) from H 

Figure 1. Dependence of K=E(du) from H. 
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K is the expectation of portfolio monthly return; H is the sensitivity of an investment rule to the factor; the 
factor is 3-month Treasury Bill rate measured in percentage terms. 
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P is the long-run covariance of the portfolio monthly return and the factor; H is the sensitivity of an 
investment rule to the factor; the factor is 3-month Treasury Bill rate measured in percentage terms. 
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Figure 3: Dependence of i? = limf^oo(l/^) Var(?i) from H 

Figure 3. Dependence of R=limt^oo(l/t)Var(u) from H. 
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R is the long-run growth rate in variance of the portfolio value; H is the sensitivity of an investment rule to 
the factor; the factor is 3-month Treasury Bill rate measured in percentage terms. 
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Figure 4: Dependence of optimal {h, i?)-strategy from risk-sensitivity parameter 

e 

Figure 4. Dependence of optimal (ii,H)-strategy from risk- 
sensitivity parameter 9. 
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r=0; H is the sensitivity of an investment rule to the factor; h is the investment amount if the deviation of 
the factor from its average value is 0; ; the factor is 3-month Treasury Bill rate measured in percentage 
terms. 
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Figure 5: Dependence of optimal {h, -ff)-strategy from risk-sensitivity parameter 
and factor-sensitivity parameter T 

Figure 5: Dependence of optimal (h,H)-strategy from risk- 
sensitivity parameter and factor-sensitivity parameter F. 
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5 Conclusion 

This paper derives explicit formulas for asymptotic joint moments of the loga- 
rithm of investor's wealth and factors, allowing for a fast maximization of linear 
combinations of the moments. An application of this method to the real data 
shows that the optimal strategy is very sensitive to the factor-sensitivity param- 
eter r that measures dependence of the investment criterion on the covariance 
of wealth with factors. 



A Proof of Theorem 1 



First, we are going to compute Eu: 



du = d{\nU) = — - -i — f 



= {h + HX)'{{a + AX)dt + Y.dW) 



(17) 




Integrating (18) over time gives 



Eu = Kt, 



(19) 



where 




14 



Our second goal is to compute E{uX). 

d{uX) = Xdu + udX + dudX 

= X[{h + HX)'{{a + AX)dt + Y.dW) -\{h + HX)'Y,Y,'{h + HX)dt\ 
+ u{BXdt + MW) 

+ [{h + HX)'{{a + AX)dt + Y.dW) -\{h + HX)'Y,Y,'{h + HX)dt]{BXdt + MW) 

(21) 

Using E{X) = E{X^) = 0, we get 

Ed{uX) = dt{AA'h + AH'a-AH'ET,'h + BE{uX)+Ai:'h) (22) 

Denoting E{uX) as P{t), we note that it satisfies the following differential 
equation 

dP 

^=BP + AA'h + AH' a - Aif'ES'/i + AE'/i. (23) 
Since B is stable, 

P{t) B-\A{H'T,T,'h - A'h - H'a) - AS'/i] as i ^ oo. (24) 

Thus, we computed P =: limt^oo P{t), which is the asymptotic covariance of 
the portfolio and factors. 
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The next step is to compute E{uXX'). 

d{uXX') = duXX' + u{dX)X' + uXdX' + dudXX' + {du)XdX' + udXdX' 

= [{h + HX)'{{a + AX)dt + T,dW) -\{h + HXyT.T.\h + HX)dt\XX' 
+ u{BXdt + MW)X' 
+ uX{BXdt + AdWy 

+ [{h + HX)'{{a + AX)dt + Y.dW) - ^{h + HX)'Y.Y.'{h + HX)dt]{BXdt + AdW)X' 
+ [{h + HX)'{{a + AX)dt + Y.dW) -\{h + HX)'Y,Y,'{h + HX)dt]X{BXdt + KdW)' 



u{BXdt + MW){BXdt + MW)' 



(25) 



Ed{uXX') = dt[h'aA + E{tr{AXX' H')XX') 

- ^/I'SS'/iA - ^E{ti{E'HXX'H'T,)XX') (26) 
+ BE{uXX') + E{uXX')B' + 2AI]'ffA + E{u)KK'] 
Asymptotically, E{uXX') Rt+S, and the coefficients R, S can be found from 
the equation: 

R = h'aA + E{tT{AXX'H')XX') - i/i'SS'/iA - ^E{tT{J:'HXX'H'i:)XX') 
+ B{Rt + S) + {Rt + S)B' + 2AE'J? A + KAA't, 

(27) 
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which reduces to two equations on R and S, 

= BR + RB' + KAA', (28) 
R = h'aA + E{tv{AXX'H')XX') - ^/I'SE'/iA - ^E{tT{J:'HXX'H'i:)XX') 
+ BS + SB' + 2AT,'HA. 

(29) 

Since A satisfies BA + AB' + AA' = 0, (28) impUes 

R = KA (30) 
and, therefore, (29) can be rewritten as 

KA = h'aA + E{tY{AXX'H')XX') - i/i'SS'/iA - ^E{ti{J:' HXX'H'J:)XX') 
+ BS + SB' + 2AE'HA. 

(31) 

Using the identity 

E{xiXjXkXi) = AijAki + AihAji + AuAjk, (32) 

we can further reduce (31) to 

BS + SB' = -2AH'AA + AH'Y^Y^'HA - 2AS'iJ A. (33) 

For later use, we also need to compute E{uX'QX) where Q is an arbitrary 
matrix. Because of the identity 

E{uX'QX) = ti{E{uXX')Q), (34) 

we have 

E{uX'QX) ~ tv{RQ)t + tr{SQ) when t oo (35) 
17 



The next step is to compute E{u'^). 

d{u^) = 2u{du) + {duf 

= 2u[{h + HX)'{{a + AX)dt + Y.dW) -\{h + HX)'Y,Y,'{h + HX)dt] 
+ {h + HX)'Y.Y.'{h + HX)dt 

(36) 

Ed{u^) = 2dt[h'aEu + E(uX')H'a + h' AE{uX) + E{uX' H' AX) 

- ^{h'Y.Y,'hE{u) + h'T,J:'HE{uX) + E{uX')H'Y,Y,'h + E{uX' H'T^Y,' HX)) 
+ E{{h + HXy^i:'{h + HX))] 

= 2dt[h'aKt + P'H'a + h'AP + tr{RH'A)t + tv{SH'A) 

- ^{h'T.Y.'hKt + h'ET,'HP + P'H'Y,Y,'h + iY{RH'Y.Y.' H)t + tr(S'ff'SS'ff)) 
+ ^E{{h + HX)'Y.i:'{h + HX))] 

(37) 

From this we have, 

Eu^ = const + t^lh'aK + ti{RH'A) - ]^{h'Y.Y.'hK + tv{RH'Y.Y,' H))\ 

+ t[2{P'H'a + h'AP + tY{SH'A)) - {h'T,T,'HP + P'H'T,T,'h + tv{SH'T,T,' H)) 

+ E{{h + HX)'Y.Y.'{h + HX))]. 

(38) 
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Finally, 

Var{u) = E{u^) - E{uf 

= const + 2{P'H'a + h'AP + tr{SH'A))t 

(39) 

- {h'HH'HP + P'H'^E'h + tv{SH'EIl'H))t 
+ (/I'SS'/i + tv{T,'HAH'E)t 
Using the definition of P this can be manipulated into 

Var{u) = - h'A - a'H)B-^K + h'T,){{h'T,T,' H - h'A - a'H)B-^K + h'T)' 

+ tr{2SH'A + (A - S)H'T,Y.'H)] + const 

(40) 
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